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Abstract 

This is an addendum to the paper “Some models of spin coherence and decoherence in 
storage rings” by one of the authors [T] in which spin diffusion in simple electron storage 
rings is studied. In particular, we illustrate in a compact way, a key implication in the 
Epilogue of [T], namely that the exact formalism of [T] delivers a rate of depolarisation 
which can differ from that obtained by the conventional treatments of spin diffusion which 
rely on the use of the derivative dnjdr] [21 [3l [!]■ As a vehicle we consider a ring with a 
Siberian Snake and electron polarisation in the plane of the ring (Machine II in [T] ). For 
this simple setup with its one-dimensional spin motion, we avoid having to deal directly 
with the Bloch equation HIE] for the polarisation density. 

Our treatment, which is deliberately pedagogical, shows that the use of dh/dr} provides 
a very good approximation to the rate of spin depolarisation in the model considered. But 
it then shows that the exact rate of depolarisation can be obtained by replacing dnjdr] by 
another derivative as suggested in the Epilogue of [T], while giving a heuristic justification 
for the new derivative. 
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Introductions 


The paper “Some models of spin coherence and decoherence in storage rings” by one of the au¬ 
thors [1] provides an introduction to the use of spin-polarisation transport equations of Fokker- 
Planck and Liouville type in electron storage rings. The work goes into great and careful detail 
and it uses two exactly solvable but simple model conhgurations called Machine I and Machine 
II as vehicles. These are perfectly aligned flat rings and only the effects of synchrotron radia¬ 
tion and of longitudinal oscillations on electron spins lying in the machine plane are considered. 
Moreover the orbit and spin dynamics are approximated by smoothed equations of motion with 
azimuth-independent parameters given by averages around the ring. Machine II contains a 
point-like Siberian Snake in addition [9]. As is often the case with “toy” models, the exact 
results for the spin distributions for Machines I and II can be expected to deliver useful insights 
and expose essential features. This is indeed the case and in this paper we use Machine II to: 


• cover a topic mentioned, but not treated in detail in [I], namely the contrast between the 
exact calculations of the rate of depolarisation in [1] and the conventional approach to 
estimating depolarisation due to Derbenev, Kondratenko and Mane (DKM) [21 [3] . 

• in the process: 

— give, and comment on, the form of the derivative dfi/drj in the DKM formula for 
this case, 

— show how to improve the DKM calculation in a heuristically obvious way, while 
conhrming a claim in the Epilogue of [1] and thereby calibrating the DKM approach 
in this case. 

• discuss approximations to dn/drj and misunderstandings in the literature. 

• suggest further avenues of investigation. 

For both model machines, polarisation build-up due to the Sokolov-Ternov effect is 

ignored. In any case the Sokolov-Ternov mechanism is ineffective in Machine II as explained 
in Section 3.2. 

In Section 2 we outline the models. In Section 3 we present the exact solution for the rate of 
depolarisation of Machine II and the rate from the DKM theory. Then, after comparing the two 
results we show how to modify the DKM calculation to get agreement with the exact solution. 
Finally, in Section 4 we make further relevant comments and discuss other approximations and 
their utility. Machine II is the simplest non-trivial exactly solvable model that we are aware 
of. 

^This paper was first prepared in 2000 following suggestions at that time that a Siberian Snake would help to 
preserve the polarisation in a putative high energy electron ring [7] . It was then put aside when the proponents 
sensibly dropped the idea. Now, after recent work [8] , it is appropriate to make an updated version of this paper 
available. 
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2 The models 


2.1 Machine I 

According to (2.5) in [1], for Machine I the Langevin eqnations of motion for the orbit and the 
spin-expectation value (the “spin”) are 

a'(s) \ / 0 -K 0\ / a{s) \ / ^ \ 

V(s) = -2as/L 0 ■ r]{s) + ■ C(s) , (2.1) 

V'Xs) / \ 0 27rz>/L 0 / \ ^/;(s) / V 0 / 

where 

• s denotes the distance around the ring (the “azimuth”), 

• L is the length of the ring, 

• cr is the distance to the centre of the bunch, 

• rj is the fractional energy deviation, 

• "^(s) is the angular position of the spin in the machine plane w.r.t. a horizontal direction 
processing according to the T-BMT equation m at the rate 27r{{g — 2 )/ 2 } 7 o/L on the 
design orbit, where {g — 2)/2 is the gyromagnetic anomaly and where 70 = Eq/{ moCo'^) 
is the Lorentz factor at the design orbit energy Eq, whereby Cq is the vacuum velocity of 
light and mo is the rest mass of the electron, 

• as is the 1-turn synchrotron damping decrement, 

• K is the compaction factor, 

• Qs = 2,71 ■ Qs/L where Qs is the synchrotron tune for undamped motion. 

• ^ = {{ 9 - 2)/2}7o, 

• a; is s-independent and is the 1-turn averaged stochastic kick strength which is expressed 
in terms of an s-independent curvature 


where in the classical electron radius, Ag is the reduced Compton wavelength of the 
electron and where is determined by the local curvature K^{s) of the design orbit in 
the horizontal plane of the original ring via the 1-turn average in the second equality, [§ 

• C simulates the noise kicks due to the synchrotron radiation, via a Gaussian white noise 
process, i.e. with stochastic averages < C('Si) • C('52) >= ^(si — S 2 ), < C('5) >= 0, where 

6 denotes the Dirac delta function. 

^In [T] w was written in Gaussian units and in terms of a product of factors Ci and C 2 - Here we work in SI 
units and express the product in terms of more convenient factors. 


U) = 


55 ,5 55 

TeKloK^ : = 


24^3 


24x/3 


TeKll- 


ds \K^ 
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More details can be found in [T] but the content of fl 2 .ip is clear: the motion of a and rj is that 
of a damped harmonic oscillator subject to the noise ( in the variable 77 . Moreover spin 
is a “passenger”, and w.r.t. the rotating reference direction it processes at the rate {2nu/L)ri. 
Then the noise acting on r] feeds through onto the spin to cause the spin diffusion which we 
wish to study. Equation (12.Ih is the cleanest way to package that message. 

In terms of the parameters 

a = -K, b = nl/K = {A'lT^QD/iKL'^) , c = -2as/L , (2.2) 


the asymptotic (s —)■ + 00 ) r.m.s widths of the distributions of a and r] are given by 


a 


2 

(7 


uia 



CO 

Yc 


cuL 

ias 



(2,3) 


We also need d = 27ri'/L. For Machine I the 1-turn periodic solution to the T-BMT equation 
on the design orbit, ho, is perpendicular to the machine plane and h is the design-orbit spin 
tune, z/q, namely the number of spin precessions around ho per turn around the ring for a 
particle on the design orbit. 

For the HERA electron ring [12] running at about 27 GeV the values of the parameters are 


Qs ^ 6.0-10-2, a, f^3.2-10-^ k 6.9 • 10-^ 
u ~ 2.0 ■ 10-^2 LRi6.3-10^m, d ^ Q.2 ■ , 

and we adopt these parameters to illustrate our arguments. Then 

a ^ -6.9 ■ 10-^ b^5.2- 10“® m- 2 , ab ^ -3.6 ■ 10"^ m- 2 , c ^ -1.0 ■ 10"® m'^ 

while 


ct2 ^ 1.0 ■ 10-® , al ^ 1.3 ■ 10"^ m2 , (2.4) 

For Machine I, the corresponding depends on the initial distribution of 7 /^ and an inter¬ 
esting point there is that the distribution of t/z reaches a stationary form after a few damping 
times. Details can be found in [D ca Ej. This is in contrast to earlier calculations which 
suggest that the 'ip distribution spreads out with a width proportional to y/s [I5] . 

However, when viewed from the hxed machine coordinates the spin distribution is not sta¬ 
tionary. On the contrary it is rotating continuously in the machine plane at the rate 2'kvq/L 
and cannot be periodic if toq ^ integer. Such a distribution cannot be handled using DKM 
methods for reasons to be explained in the Commentary. 

2.2 Machine II 

For Machine H a single point-like Siberian Snake is included at s = 0. This rotates a spin by 
the angle vr around the radial direction independently of a and p, so that it is a “Type 2” snake 
eg, In Machine I, ho, is perpendicular to the machine plane. But for Machine H ho is in the 
horizontal plane and is given by (3.3) in [1] H: 

ho(s) = cos( 5 ( 0 ( 5 ) lei-h sin( 5 ( 0 ( 5 ) )62 , (2.5) 


^For clarity we omit the subscript II used in [T]. 
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where ei and 62 are unit vectors transverse and parallel to the design orbit respectively and 

= d(^s-L/2-Lg{s/L)^ , (2.6) 

where the “stairway” function Q is defined by: 

g{s/L) = N a N < s/L < {N + 1) , (2.7) 

for which N is an integer. In the range 0 < s < L one has 

geis) = d(^s - L/2^ . ( 2 . 8 ) 

For Machine II it is this horizontal ho(s) which is taken as the reference direction so that 
for Machine II is the angle between a spin and ho(s). The conclusions of this paper are the 
same for a snake with a longitudinal rotation axis (“Type 1”), provided, for example, in the 
real world, by a solenoid. In that case the spin motion on the design orbit is illustrated in 
Figure 8 in [16]. The reader can easily make a corresponding sketch for our Type 2 snake. 

Owing to the presence of the snake, z/q is not u but 1/2 [Il[l 6 |. Thus the condition spin-orbit 
resonance occurs around Qs = 1/2. Near resonance, spin diffusion effects can be particularly 
strong. 0 


3 The polarisation evolution for Machine II 


3.1 Using exact distribution functions 

We now look again at the exact evaluation of the distribution functions for Process 3 of Machine 
11 [D Section 3.5]. For this process the orbital phase-space distribution is in its stationary state 
and the spins are initially all parallel to ho. Then by (3.58) and (3.62) in [T], and after transients 
have died away in the first few damping times the polarisation of the whole beam evolves like 


Pm{s)\ \ oc exp 


9uis) \ 
2 ; 



(3.1) 


where the function 


gi 4 ^ is 1 -turn periodic and gi^, which is positive, is given by 

2 Oil / \ 

^15 = - , j ' ( 2A sinh(cL/ 2 ) - csin(AL) ) , 

a b L X \ ) 


(3.2) 


with A = \/—ab — jA and gu = l/{cosh(cL/ 2 ) -|- cos(AL)}. 

Whereas Aq ^ \/—ab = ilg is proportional to the synchrotron tune in the absence of damp¬ 
ing, A = ^J—ab — c^/4 is the corresponding tune in the presence of damping. Since « —ah, 
A Aq. 

At large times the polarisation vanishes: 

ll^to/(+oo)|| = 0. (3.3) 

^The resonance condition should be more precisely expressed in terms of the amplitude-dependent spin tune 
[min]. But for typical electron/positron rings the amplitude-dependent spin tune differs only insignificantly 
from vq. 
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The depolarisation rate w.r.t. distance is: 


(3.4) 


1 _ 9i5 

T~spin 2 

For the HERA electron ring parameters listed above one gets 

Tspin ~ 7.6-lO^m, (3.5) 

which corresponds to about 12000 turns, i.e. about 250 milliseconds. 

For the HERA parameters, and when transients have died away after about 1000 turns 
(s ~ 6.3 • 10® m), the variation of gui^s) over a turn causes a variation of 5 'i 4 (s) + sgi^ in the 
exponent of fl3.ip of about 15 percent. If Qg were close to 1/2 so that XL were close to tt, 
5^15 would, because of its factor gn, become very large and Tgpin would be very small. This is 
exactly what one expects when sitting close to a spin-orbit resonance HT]. 

Of course this model only includes the effects on spin of smoothed synchrotron motion and 
radiation in the main body of the ring. No account is taken here of the detailed dependence 
of spin motion on the orbital variables in a real snake and there is no horizontal or vertical 
betatron motion. 

It can be seen from [1] that the asymptotic depolarisation rate is ( 715/2 for any initial 
distribution of spins in the machine plane. 

For comparison with the results of the next section we use the relations = —u:/2c and 
Ao = V—ab to obtain 

Apin = 4 • 9 r ■ 7 - r/o! I - am ■ f^A sinh(cL/2) - csin(AL)^ . (3.6) 

^ Aq 2c XL {cosh(cL/2)-|-cos(AL)| \ / 

3.2 Using the Derbenev—Kondratenko—Mane approach 

The conventional way to calculate the rate of depolarisation is to use the spin diffusion term 
in the Derbenev-Kondratenko-Mane formula PE] for the equilibrium polarisation. In this 
approach it is assumed that at orbital equilibrium, the combined effect of the depolarisation 
and the S-T mechanism is to cause the polarisation at a point in phase space to be aligned along 
the vector h of the invariant spin field [UllIH]. This is a special solution to the T-BMT equation 
along the trajectory (cr(s), 77 ( 5 )) satisfying the periodicity condition h{a, g; s) = h{a, g] s + L). 
In general h is a function of all six phase-space coordinates but in our models only a and g come 
into play. In the absence of radiation, and spin-orbit equilibrium, the polarisation at a point 
in phase space, which we call the local polarisation is indeed aligned along fi [19]. Then, since 
the characteristic times for the action of synchrotron radiation, namely the damping time and 
the polarisation and depolarisation times, are very large compared to the characteristic times 
for the orbital and spin-precession dynamics, the above assumption about the direction of the 
local polarisation for electrons with radiation is reasonable since the orbital and spin dynamics 
still dominate on short time scales, at least away from spin-orbit resonances. We return to this 
in Section 3.4. The DKM approach also assumes that the value of the local polarisation and 
its rate of change are independent of the position in phase space. These two assumptions are 
motivated by recognition that the stochastic photon emission and damping cause electrons to 
continually diffuse through phase space and thereby, in the end, have effectively interchangeable 
histories. 
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According to the DKM approach, with a stable phase-space distribution, the local polari¬ 
sation settles to the equilibrium value 


P. 


f ds ■ (h 


dh ' 
dr ;' 


dkm 


5^3 


f ds ( |A:(s)| 3(1 - |(h -sf + % 


(3.7) 


where s = direction of the particle motion, b = [s x s)/|s|, K{s) is the curvature (in the 
horizontal or vertical planes) and < >s denotes an average over phase space at azimuth s. 
The unit vector b is the magnetic held direction if the electric held vanishes and the motion is 
perpendicular to the magnetic held. The polarisation of the beam as a whole is 

^dkm(s) = Pdkm {fl)s (3.8) 

In the DKM formula, the depolarisation rate w.r.t. time is 



dep 





1 

L 


ds 





(3.9) 


where, by the very nature of the DKM approach, it is assumed that all transients have died 
away. By (13.8p this is also the depolarisation rate of the whole beam. For a ring with constant 
curvature fl3.9p becomes 



dep 


55 

48\/3 


J^3 

reXeCol^^ 


ds 


so that the depolarisation rate w.r.t. distance is 



(Co ‘ = 



(3.10) 


To evaluate this we need to know h at each azimuth and at each point in phase space. From 
(B.l) in Appendix B in [T], in the arc the T-BMT equation reads as 


dh 

ds 


- 0 - V ^ - b ^ + d{l + 77)63 X h . 
da dr] 


(3.11) 


This is in fact just the partial differential form of the T-BMT equation 


dh da dh dr] dh 
ds ds da ds dr] 


d{l + 77)63 X h 


(3.12) 


for which the derivatives da/ds and dr]/ds have been obtained from fl2.1l) by switching off the 
radiation. 

By choosing the ansatz 


h = cos(/)ei + sin(/)e 2 


^For clarity we still omit the subscript II used in [T]. 


(3.13) 









fl3.1ip becomes 


(3,14) 


ds 


df.df 

-a r] - - ba — + dr] + d. 

oa or] 


Then by taking the snake into account and enforcing 1-turn periodicity we obtain 


f{a ,)); s) = ge(s) + agi,{s) + ggxis) , 


(3.16) 


with 


9i9{s) = 


d b 


Ag 1 T cos(AqT) 


cos( Ao[s — L — LQ{s/L)] 


+ cos( Ao[s — LQ{s/L)] j — cos(Ao-h) — 1 


92o{s) = 


d 


Aq 1 T cos(AoT) 


sin( Xo[s — L — LQ{s/L)]] + sinf Aq[ s — LQ (s/ L)] 


(3.16) 


where both are independent of a and rj. 

The resulting h is a solution of the T-BMT equation along the trajectory {a{s),ri{s)) and 
like fiQ it lies in the machine plane. As required, it is 1 -turn periodic in s for all a and rj and 
reduces to ho at a = 0, = 0. If n had a component perpendicular to the machine plane it 

would, at simplest, be 2 -turn periodic. 

One also sees that a singularity in h occurs if the fractional part of the orbital tune Qg = 
(Ao ■ L)/{27i) equals 1 / 2 , i.e. if one is at a spin-orbit resonance dehned in terms of z/q (see 
Footnote 3). In that case h(cr, s) is not unique and a different formulation is needed. 

Obviously, [dh/dr])^ = which is independent of a and rj. Then to obtain the rate of 
depolarisation within the DKM framework we just need to evaluate 


(co T-dep) - ^ 


glois') ds'. 


In the range 0 < s < L we have 


92 o { s ) = 


d 1 


Ag 1 T cos(AoT) 


sini Ao(s — L) ) + sin(Aos) 


so that we hnd 


(co 


V-l 


OJ 


d^ 


Ag^ 2 Ag T {1 + cos(AgT)} 


AgL — sin(AgL) ) . 


(3.17) 


(3.18) 


(3.19) 


This scales like 7 / at hxed Ag. Thus the depolarisation due to noisy damped synchrotron 
motion becomes very strong at high energy. 

For Machine II, the vector h is in the machine plane so that b ■ n is zero and the 
S-T mechanism is inoperative |lj. However, in this case the vector dh/drj acquires a vertical 
component. This leads to the phenomenon of kinetic polarisation embodied in the term linear 
in dn/drj in the numerator in fl3.7p [201 [T 6 ] . We will not pursue that here. The formalisms 
in PE] leading to fl3.7p and fl3.9p are semiclassical and the concept of depolarisation does 
not immediately appear in those papers. Nevertheless, the expression for in 03.91) can be 
obtained from classical notions as in [2T1|22|. This fact will be useful in later discussions. 
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3.3 Comparison 

Since c L k, —6.3 ■ 10“^, and A ~ Aq it is immediately clear that the result for (cq T^ep)”^ in 
(I3.19P is very close to in (I3.6p and in fact the relative difference between the two rates is 
^ 4 - 10 -®. 

Thus the DKM estimate, that used is all analytical calculations of the depolarisation rate, 
is very accurate and perfectly adequate for Machine II with the parameters used here. In fact 
it is expected to be very accurate in general. But it is not exact. This is easy to understand: 
the DKM result is based on the impressive insight that the equilibrium polarisation should be 
parallel to n at each point in phase space but as we will see in the next section the polarisation 
is not exactly parallel to n. 

If we set Qs close to 1/2 so that we are close to spin orbit resonance [§ the DKM estimate 
deviates signihcantly from the exact result. For example for Qs = 0.4998 the relative differ¬ 
ence is about 65%. We expect this to be the case in general for realistic rings too. But the 
depolarisation rate would then be so large that the results would be of no practical interest. 


3.4 A modified DKM calculation 


Given the (admittedly small) differences between the two values for the depolarisation rate we 
are motivated to explain them by exploiting our exact expressions in [1] for the asymptotic 
spin distributions and thereby examine the assumptions in Section 3.2 underlying the DKM 
approach. 

Thus, as shown by (2.120) in [T], the unit vector P^^^{a,ri;s) describing the polarisation 
direction at each point in phase space does not(!) obey the T-BMT equation. In particu¬ 
lar, damping and diffusion effects must be included. Then at orbital equilibrium and after 
transforming (3.33) in [1] into the machine frame we hnd 


dP: 


dir 


ds 


= —a T] 


dP, 


dir 


dP, 


da 


a 


dir 


Orj 


iljj X 


pw 

dir 


+ 0 7 ] 


dP, 


dir 


dr] 


(3.20) 


It was shown in [T] that after transients have died away in the hrst few damping times, 
P^^^{a,r];s) becomes 1-turn periodic is s although the local polarisation itself then decreases 
smoothly to zero. From now on we will denote the periodic P^^^{a,r]] s) by the unit vec¬ 
tor p{a, 7]] s) and we will now show how to obtain it in a way paralleling the construction of 
n{a,r]-, s). According to fl3.20p . in the arc p fulhlls 


dp 

ds 



UjJ UjJ . , „ 


(3.21) 


Then, writing 


p = cos(/)ei -F sin(/)e 2 


(3.22) 


we have, in contrast to (13. Mh . 


ds 


df , 

-a h ^ - 0 0- 


df df 1 1 

^ + cp — + dp + d. 
dp dp 


(3.23) 


®Of course one never tries to run a storage ring in that way — the r.f. cavity voltage would be enormous 
and the smoothed equations of motion used in the model would no longer be reasonable. 
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So there is an extra term depending on the damping rate. Taking into account the action of 
the snake, the 1 -turn periodic solution for / is 

/(a, ? 7 ; s) = gais) + < 7 ^ 19 ( 5 ) + vhois) (3.24) 

where in the range 0 < s < L 


^19(5) = 


d 


2A a 


i 9i{s) gii exp(-c L/2) + i gi{s - L) gu exp(c L/2) - 2 A 


(3.25) 


^2o(s) = 


d 


2 A 


i g 2 {s) gii exp(-c L/2) + i g 2 {s - L) gu exp(c L/2) 


(3.26) 


with 


gi{s) = i exp(c s/2) [c sin(A s) — 2 A cos(A s)] , 

g 2 {s) = 2i sin(A s) exp(c s/2) , (3.27) 


so that fl3.22p and (I3.24p are equivalent to (3.72) in [T]. Thus p{a, Tj] s) 7 ^ h(cr, 77 ; s) owing to 
the dependence of p on c but for c —)■ 0 we see that ^ 19 ( 5 ) —)■ giQ^s), P 2 o(s) —)■ ^(s) and p ^ n 
as one would expect. Moreover, the expression for T~p]^ in (13.bh then reduces to the expression 
for (co fl3.16p . The function g 2 o{s) can be written in the form 


g2o{s) — - gu 


sin(A s) exp(c (s — L)/2) + sin(A (s — L)) exp(c s/2) 


(3.28) 


and apart from the exponential factors containing c, it is reminiscent of g 2 Q in fl3.18p . 

We now turn to the local polarisation s)\\. In the DKM picture it is assumed 

that at orbital equilibrium | |T))('))(cr, p; s)|| and its rate of change are independent of {cr,p). As 
seen in (3.71) in [1], this is indeed the case in our model after transients have died away. This 
emerges naturally - there has been no need to appeal to heuristic arguments. In particular, 
at orbital equilibrium the long-term s dependence of ||T)of(s)|| reflects the s dependence of 
s)\\. However the direction of the local polarisation is p, not h. Processes 4 and 5 in 
[ 1 ] reflect on this too. 

As we are seeing, instead of relying on intuition, for Machine II one can hnd the exact 
result for the asymptotic depolarisation rate by examining the development of the distribution 
functions. In the general case, such an analytical treatment is not available. But to stay close 
to the philosophy of using distribution functions one can examine the properties of solutions of 
the Bloch equation for the polarisation density [B O | 6 ] and then, as suggested in the Epilogue 
in [1], (dh/dp)'^ should be replaced by {dp/dp)"^ in the expression for (cq T^ep)”^ of DKM 
picture. It is then no surprise to hnd that 


L 

(cor.^p)"^ = ^ j g^i^') ds'. 
0 


(3.29) 


gives precisely the in fl3.6p ! 
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This is also exactly what one expects given the physical picture in [2T1122] used to arrive at 
the use of {dn/drjY if the polarisation were parallel to n. There, a simple geometrical construc¬ 
tion was used to estimate the rate of change of the average of spin projections along the local 
polarisation direction, which was taken to be n. For a stored non-radiating beam in a station¬ 
ary spin-orbit state, i.e. a non-radiating beam for which the phase-space distribution and the 
polarisation distribution repeat from turn to turn after a chosen starting s, the polarisation 
is indeed parallel to n at each point in phase space im. But in the presence of radiation the 
direction of the polarisation is p and if p had been used in [211 122 ] instead, then {dp/dpy would 
have been needed as we have already discovered by appealing to the Epilogue in [1]. Thus the 
exact agreement between (I3.29p and (I3.6p supports the construction used in [2ll |22] once h has 
been replaced by p. The necessity of using p instead of h emphasises that the DKM approach 
is essentially perturbative beginning with n. 

Of course, the “improvement” embodied in fl3.29p is so small that it is of no practical 
significance. Nevertheless our discussion demonstrates how our toy model can be used to 
expose details of spin dynamics which would not otherwise be accessible. 

As seen in [1] , at the phase-space point where a = and rj = arj and for Qg well away from 
1/2, the angle between n and p at the snake is given approximately by 

d ■ c - ar, \o ■ L — sin(Ao ■ L) 

2 • Aq 1 -1- cos(Ao • L) 

For HERA parameters this is about 0.04 milliradians whereas the angle between n and 
ho is of the order of 200 milliradians. The angle between h and p is proportional to c so 
that, as one might expect, the angle becomes larger as stronger damping causes the equation 
of motion for p to deviate more from the T-BMT equation. Moreover, by comparing the 
denominators in and (I3.16p we see that p differs strongly from h near Qg = 1/2, thereby 
confirming our suspicion in Section 3.2 that near spin-orbit resonances, care is needed with the 
argument based on time scales. This, then, is the origin of the failure of the near-resonance 
DKM estimate for the depolarisation rate mentioned in Section 3.3. Note that as seen in [1] 
the characteristic time for the decay of transients remains the damping time, even close to the 
resonance. In [23] it is suggested that very close to spin-orbit resonances, extra terms containing 
delta functions should be added to the DKM value (I3.9p for the rate of depolarisation as a result 
of so-called uncorrelated resonance crossing. These would be due to fluctuations in the rate 
of spin precession around ho caused by energy fluctuations. See [8]. However, for Machine H 
spins precess only around the vertical in the arc, not around ho which is horizontal in the arc. 
So this effect will not be seen here. 


4 Commentary 

Some further remarks are now in order —some obvious, some less so. 

(1) It is clear from the discussion above that ho and h are two different quantities which only 
coincide on the closed design orbit. Indeed, as we have just mentioned, the angle between 
them is of the order of 200 milliradians for our parameters. However, they have often been 
confused in the literature. For example dn/drj was originally written as 'y{dn/d'y) [JT] |2] 
and that led some to calculate 'yo{dno/d'yo). See, for example [16]. For general problems 
where all six phase-space coordinates must be included this can give completely misleading 
results. In particular there would be no resonant increase of the depolarisation rate when 
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an orbital tune were close to z/q. Nevertheless, if, as in our models, horizontal and vertical 
betatron motion are being neglected, it can happen that (7o<9ho/(97o)^ provides a useful 
initial approximation to [dfi/dr])'^. That is the case with Machine II as we now show. 

Using fl2.5p and fl2.8p and the definition of d and u 


= d^ {s- L/2f = {v {2 -k s/L - Tr)}^ . (4.1) 

d7o d7o 


This is to be compared with 






2 

{1 + cos(AoT)} 


• sin^Ao(s — L/2) 


( 4 . 2 ) 


The expression {v (27r s/L — vr)}^ in fl4.ip is used in [20] and its origin is clear from the 
absence of the resonance factor 2/{l + cos(AoT)} which takes the value ~ 1.036 instead 
of 1 for our value: Qs = 6.0 • 10“^. By replacing {dn/drj)‘^ in (I3.19p by (7o<9ho/(97o)^ and 
using the result 


L 



0 


one obtains a value about 3% lower than from fl3.19p so that in this case (7o<9ho/(97o)^ 
provides an adequate approximation. But of course it would become a bad approximation 
for a large Qs- Note that the approximation works for Machine II because the snake 
ensures that the spin tune, vq, is far from Qg. In other situations one should never rely 
on this approximation. The use of (7o<9ho/57o)^ delivers the correct result for Machine II 
at Qs = 0, i.e. when the energy is the same from turn to turn. 

(2) The vector h is a s-periodic solution to the partial differential equation (13.lip . Whereas 
ho can be obtained as the unit real eigenvector of the 1-turn spin map on the closed 
orbit [16], h is not the eigenvector of the 1-turn spin map beginning at some a, ?7 and 
s unless Qs is an integer. At each chosen hxed s, each of the three components of h 
lies on a component-specific closed curve in the (a, rf) plane corresponding to the closed 
ellipse in the (a, rf) plane mapped out by a non-radiating particle of fixed amplitude. But 
if a particular (a, rj) pair are chosen at some initial s and the corresponding h((T, rj] s) 
is transported according to fl3.1ip for one turn, this h does not in general return to its 
starting value so that it is not a “closed spin solution”. 

(3) Some authors still use the symbol n when they actually mean ho! The tendency to create 
confusion seems to be deep rooted. 

(4) The vector dfi/drj (which is still often written as 'jdn/d'y) is sometimes called the “spin 
chromaticity”. We prefer the terms “spin-orbit coupling function” or “spin-field deriva¬ 
tive” so that “spin chromaticity” can be reserved for the rate of change of a amplitude- 
dependent spin tune w.r.t. a fractional energy change. In any case in the full theory, the 
DKM formula fl3.7p must be modified to include (usually) relatively small terms involving 
derivatives of h w.r.t. the two transverse canonical momenta [HI 121] and for such terms 
the name “chromaticity” is clearly unsuitable. 
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(5) If IcT^fig + rig 2 o\ « 1, i.e. if the angle between h and ho is small, then h in fl3.13p and 
(13.1 Sp can be approximated by 


where 


h = ho + 0-5(19(5) + r]g2o{s) m 


rh = - sin ( g^^s) ) ei + cos ( g^^s) ) 62 


(4.4) 

(4.5) 


is a unit vector in the machine plane perpendicular to h. This is the “SLIM approxima¬ 
tion” [25] whereby h deviates from hg by a function linear in the phase-space coordinates 
[Ij. In this approximation we again hnd {dn/drjY = 5'2o Machine II the 

SLIM approximation actually delivers the correct result for [dn/dr])'^. Of course, a direct 
application of the eigenvectors and matrices of the SLIM formalism delivers (14.4p too. 

(6) The reason that the DKM approach cannot be used on the spin distribution of Machine 
I is as follows. The DKM formalism assumes that transients have died away and that the 
polarisation is locally parallel to h. But in Machine I the magnetic held is perpendicular 
to the machine plane everywhere so that tiq and n are also perpendicular to the machine 
plane. Then since the spin distribution of Machine I is set up in the machine plane the 
spins are perpendicular to h and the distribution processes at the rate z/q around Uq. In 
Machine II there is a magnetic held in the machine plane at the snake and uq and h are 
in the machine plane together with the spins themselves. In this case the polarisation can 
settle down to be almost parallel to h. 

(7) We have seen how the holds h((T, rj] s) and p{a, rj] s) for Machine II are found by enforc¬ 
ing 1-turn periodicity while solving the partial diherential equations (PDE) (I3.1ip and 
(I3.2ip respectively. For real rings, the most reliable method for hnding n is stroboscopic 
averaging [26] for a point in phase space. That involves integrating the T-BMT equation 
for spins along a particle trajectory through that point in phase space and averaging the 
spins. The 1-turn periodicity emerges automatically from the algorithm. For (13.lip the 
trajectory is, of course, dehned by the relations da/ds = ag and dg/ds = ba. But these 
just dehne the characteristic curves of the PDE, familiar from the method of characteris¬ 
tics for solving linear PDF’s [27]. Stroboscopic averaging for h relies on the fact that the 
T-BMT equation (equivalently (I3.11IB is linear. The equation of motion for p in (I3.2ip 
and for real rings is also linear. Thus, since stroboscopic averaging works well for h, it 
is tempting to use it for p. In that case one would again integrate along a characteris¬ 
tic curve. However, the dehning equations for a characteristic curve would now include 
damping. Then care would be needed. In any case stroboscopic averaging is only useful 
if the average is normalisable. 

(8) In principle the asymptotic direction of the local polarisation at points in phase space 
and its dependence on the distance to spin-orbit resonance could be discovered with a 
Monte-Carlo spin-orbit tracking simulation [28118]. This would require a large number of 
particles and corresponding computing power in order to sufficiently populate a sufficient 
number of points in phase space. In addition a precise numerical determination of h at 
each point in phase space would be needed for comparison and that would require strobo¬ 
scopic averaging [26] or the SODOM algorithm [29]. Away from spin-orbit resonances the 
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required number of particles could perhaps be reduced by just looking (say) turn-by-turn 
at the projection of a spin on the plane perpendicular to n for each particle and then 
hnding the average projection. This would represent the deviation of the direction of the 
local polarisation from h, averaged over phase space. 

(9) Although it would be satisfying to calculate the rate of depolarisation for real elec¬ 
tron/positron machines by integrating the Bloch equation for the polarisation density, 
that would be cumbersome in practice. In any case the DKM approach, as well as be¬ 
ing elegant, provides an extremely good approximation perfectly suited to the operating 
regimes of storage rings in operation up to now and it provides the only practical analyt¬ 
ical approach for real rings with arbitrary discrete magnet structures. Computer codes 
which evaluate the DKM formula numerically for real rings are listed in jl]. Unfortu¬ 
nately, beyond the hrst order approximation of SLIM, they all require large amounts of 
computing time. 

Nevertheless, the Bloch equation for the polarisation density, augmented by a Baier- 
Katkov-Strakhovenko expression for the influence of the Sokolov-Ternov effect O |6l l30] 
remains the key to a general description of polarisation dynamics in electron/positron 
rings. Note that as explained in [5l [6] and hinted at in the Epilogue in |T], with this 
combination it will be possible to arrive at the DKM estimates from hrst principles. A 
hrst approach, based on a semiclassical calculation, can be found in [23] . This programme 
will involve approximations but will enable analytical exploration of the limitations of the 
DKM formula without recourse to heuristics. See [8] for Monte-Carlo simulations of the 
depolarisation process at very high energy and for a study of whether the extra terms for 
the depolarisation rate mentioned in Section 3.4 are needed. 


Conclusion 

Analytical calculations for the polarisation in electron/positron storage rings are usually based 
on the DKM formalism. This, in turn, is based on some implicit reasonable considerations 
of the various time scales involved in spin-orbit motion. See, for example. Figure 1 in [TB] . 
Nevertheless there are open questions about the applicability of the DKM formula at very high 
energy [8]. 

It therefore seems desirable to return to hrst principles and obtain the rate of depolarisation 
from a study of the solutions of the Bloch equation for the polarisation density m- To obtain 
the equilibrium polarisation, the Sokolov-Ternov ehect must be included via the Baier-Katkov- 
Strakhovenko formalism. 

This paper has used a simple, non-trivial model of spin dihusion to show that the DKM 
expression for the depolarisation rate is only a good approximation in this case. Although the 
Bloch equation for the polarisation density was not used directly, its use was implicit and our 
results support the claim in the Epilogue of [T] that the derivative On/dr] should be replaced by 
dp/dr]. Nevertheless, further work will be needed to properly establish the range of applicability 
of the basic DKM formalism and to address the open questions mentioned in [8|. 
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